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1 Introduction 

Transport theory is a very convenient tool to study many-body nonequilibrium 
systems, nonrelativistic and relativistic as well. The kinetic equations which play 
a central role in the transport approach can be usually derived by means of simple 
heuristic arguments similar to those which were used by Ludwig Boltzmann over 
hundred years ago when he introduced his famous equation. However, such arguments 
fail when one studies a system of very complicated dynamics as the quark-gluon 
plasma governed by QCD. In spite of some efforts, the general structure of collision 
terms in the transport equations of the quark-gluon plasma remains unknown. In such 
a situation one has to refer to a formal scheme which allows to derive the transport 
equation directly from the underlying quantum field theory. The formal scheme is 
also needed to specify the limits of the kinetic approach. Indeed, the derivation shows 
the assumptions and approximations which lead to the transport theory, and hence 
the domain of its applicability can be established. 

The so-called Schwinger-Keldysh |]l[ formulation of the quantum field theory pro- 
vides a very promising basis to derive the transport equation. Kadanoff and Baym 
developed the technique for nonrelativistic quantum systems, which has been further 
generalized to relativistic ones 0, |^, |, ^ 0, ||, ^, |10|1- We mention here only the papers 
which went beyond the mean field or Vlasov approximation and provide a more or 
less systematic analysis of the coUision terms. 

The treatment of the massless fields, which are crucial for the gauge theories 
as QED or QCD, is particularly difficult when the transport equations are derived. 
Except the well known infrared divergences which plague the perturbative expansion, 
there is a specific problem of nonequilibrium massless fields. The inhomogeneities in 
the system cause the off-mass-shell propagation of particles and then the perturbative 
analysis of the collision terms appears hardly tractable. More specifically, it appears 
very difficult, if possible at all, to express the field self-energy as the transition matrix 
element squared and consequently we loose the probabilistic character of the kinetic 
theory. The problem is absent for the massive fields when the system is assumed 
homogeneous at the inverse mass or Compton scale. This is a natural assumption 
within the transport theory which anyway deals with the quantities averaged over a 
certain scale which can be identified with the Compton one. 

The problem of the massless nonequilibrium fields has not been fully recognized 
in the existing literature. One has usually assumed, explicitly or implicitly, the on- 
mass-shell propagation or equivalently the ideal homogeneity of the system. We find 
such an approach unsatisfactory and suggest a systematic solution of the problem. 
The basic idea is rather obvious. 

The fields which are massless in vacuum gain an effective mass in a medium 
due to the interaction. Therefore, the minimal scale at which the transport theory 
works is not an inverse bar mass, which is infinite for massless fields, but the inverse 
effective one. The staring point of the perturbative computation should be no longer 
free fields but the interacting ones. In physical terms, we postulate existence of the 
massive quasiparticles and look for their transport equation. 



At technical level, we begin with the lagrangian of the massless fields and make a 
formal trick which is well known in the quantum field theory at finite temperature, see, 
e.g., []12|, |13|, |14|, |1^, |1^, 0, [1^, |T^. Namely, the auxiliary mass term is added to the 



free lagrangian and then is subtracted due to a redefinition of the interaction one. As 
a result the subset of diagrams which contributes to the mass, which is determined 
in a self consistent way, is effectively resummed in the perturbative expansion. A 
somewhat similar technique was applied in the kinetic theory in [10|. 



In this paper we show how the suggested method works for the self-interacting 
scalar fields. We discuss in detail the (j)^ and (p^ models which appear to be quali- 
tatively different. We successfully derive the transport equations for the (f)'^ model 
and show why the method does not work for the 0^ case. Our discussion closely 
follows the scheme of derivation which was earlier developed for the massive fields: 
self-interacting scalar fields [0 and the spinor fields interacting with the scalar and 
vector ones |§]. 

The main steps of the derivation are the following. We define the contour Green 
function with the time arguments on the contour in a complex time plane. This 
function is a key element of the Schwinger-Keldysh approach. After discussing its 
properties and relevance for nonequilibrium systems, we write down the exact equa- 
tions of motion i.e. the Dyson-Schwinger equations. Assuming the macroscopic 
quasi-homogeneity of the system, we perform the gradient expansion and the Wigner 
transformation. Then, the pair of Dyson-Schwinger equations are converted into the 
transport and mass-shell equations both satisfied by the Wigner function. The later 
equation allows one to identify the initially introduced fictitious mass with the effec- 
tive one generated by the interaction. We further perform the perturbative analysis 
showing how the Vlasov terms and the coUisional ones emerge. Finally we define the 
distribution functions of standard probabilistic interpretation and find the transport 
equations satisfied by these functions. 

Throughout the article we use natural units where h = c = 1. The signature of 
the metric tensor is (+, — , — , — ). As far as possible we keep the convention of Bjorken 



and Drell |0| 



2 Preliminaries 

We consider the system of massless scalar fields with the lagrangian density of 
the form 

C{x) = ld'^<j){x)d,<j)ix) - ^r{x) , (1) 

2 n\ 

where n equals 3 or 4. The renormalization counterterms are omitted in the la- 
grangian. We introduce an auxiliary position dependent mass in^{x) which can be 
treated as an external field. Specifically, we redefine the lagrangian as 

Cm{x) = ]^d^<P{x)d^<P{x) - Kll{x)<P\x) + Ci{x) 



with the interaction term 

Cj{x) = +l-ml{x)<f)\x) - ^r{x) . 
2 n! 

The fields which satisfy the equation of motion 

[d' + ml{xMx) = 0, (2) 

represent free quasiparticles with mass m^. We observe that it is not a priori clear 
whether massive quasiparticles emerge due to the field self-interaction. It is even less 
clear whether the limit of free quasiparticles exist. As will be shown it is indeed the 
case for the 0^ model, but not for the 0^ one. 

We write down the energy-momentum tensor defined as 

T'^'ix) = d''(j){x)d''(j){x) - g^^Ciyx) . 

Subtracting the total derivative 

jd^d''(t)^{x)-g^'''-d''d„<P'^{x) , 

we get the energy-momentum tensor which for the free fields is of the form convenient 
for our purposes i.e. 

To'^ix) = -\<P{x) d'd^ 0(x) . (3) 

The fields are assumed here to satisfy the equation of motion (^). 

3 Green functions 

The central role in our considerations plays the contour Green function defined 

as 

iAix,y)''^{f<p{x)<Piy)) 

where (...) denotes the ensemble average at time to (usually identified with — cxo); T is 
the time ordering operation along the directed contour shown in Fig. 1. The parameter 
tmax is shifted to -|-cxd in the calculations. The time arguments are complex with an 
infinitesimal positive or negative imaginary part, which locates them on the upper or 
on the lower branch of the contour. The ordering operation is defined as 

f(j){x)(j){y) = e(xo,i/o)0(a;)0(y) + Q{yo,xo)(j){y)(j){x) , 

where 6(xo,yo) equals 1 if Xq succeeds yo on the contour, and equals when xq 
precedes yo- 

If the field is expected to develop a finite expectation value, as it happens when 
the symmetry is spontaneously broken, the contribution (0(a;))(0(?/)) is subtracted 



in the right-hand-side of the equation defining the Green function, see e.g. 0, ||. 
Then, one concentrates on the field fiuctuations around the expectation values. Since 
(0(x)) is expected to vanish in the models defined by the lagrangians (|1|) we neglect 
this contribution in the Green function definition. 

We also use four other Green functions with real time arguments: 

iA>{x,y) = {(t){x)(t){y)) , 

iA<(x,?/) = {(j){y)(f){x)) , 

tA%x,y)''^{r<P{x)<P{y)), 

tA%x,y)''^{T^<P{x)<P{y)), 
where T'^{T°') prescribes (anti-)chronological time ordering: 

T^0(x)0(y) = e(xo - yo)0(a;)0(y) + 0(yo - Xo)(l){y)(j)c{x) , 

T''(t){x)(j){y) =^ 0(|/o - xo)0(x)0(y) + e(xo - yo)0(l/)0(a;) . 
These functions are related to the contour Green functions in the following manner: 

A.^{x,y) = A{x,y) for xq, y^ from the upper branch, 

A°'{x,y) = A{x,y) for xq, yo from the lower branch, 

A-^{x,y) = A{x,y) for xq from the upper branch and 

yo from the lower one, 

A^{x,y) = A{x,y) for xq from the lower branch and 

yo from the upper one. 

One easily finds the identities which directly follow from the definitions 

A'{x, y) = 0(xo - l/o)A>(x, y) + 6(^0 - a;o)A<(x, y) , (4) 

A'^(a;, y) = 6(^/0 - 2;o)A>(x, y) + Q{xo - yo)A<ix, y) . 
One also observes that 

(iA^{x,y)^ = iA^{x,y) , 

(zA«(a;,y))^=2A^(a;,y), 

where f denotes hermitian conjugation, i.e. complex conjugation with an exchange 
of the Green function arguments. Because the fields are real, the functions iA^{x,y) 
satisfy the relation 

A>{x,y) = A<{y,x) . (5) 



It appears convenient to introduce the retarded (+) and advanced (— ) Green 
functions 

A^{x,y)''^' ±{A>{x,y) - A<{x,y))Q{±XoTyo) ■ (6) 

One immediately finds the identity 

A+(x, y) - A-{x, y) = A>(x, y) - A<(x, y) . (7) 

Let us now briefly discuss the physical interpretation of the Green functions. The 
function A'^[x,y) describes the propagation of disturbance in which a single particle 
is added to the many-particle system in space-time point y and then is removed from 
it in a space-time point x. An antiparticle disturbance is propagated backward in 
time. The meaning of A°'{x,y) is analogous but particles are propagated backward 
in time and antiparticles forward. In the zero density limit A^(x, y) coincides with 
the Feynman propagator. 

The physical meaning of functions A^{x,y) and A'^{x,y) is more transparent 
when one considers the Wigner transform defined as 

A^{X,p) = Jd\e'P^A^{X + lu,X- in) . (8) 

Then, the free-field energy-momentum tensor @ averaged over ensemble can be ex- 
pressed as 

One recognizes the standard form of the energy-momentum tensor in the kinetic the- 
ory with the function iA'^{X,p) giving the density of particles with four-momentum 
p in a space-time point X. Therefore, iA^{X,p) can be treated as a quantum ana- 
log of the classical distribution function. Indeed, the function iA^{X,p) is hermitian. 
However it is not positively definite and the probabilistic interpretation is only approx- 
imately valid. One should also observe that, in contrast to the classical distribution 
functions, zA^(X, p) can be nonzero for the off-mass-shell four-momenta. 

4 Green Function Equations of Motion 

The Dyson-Schwinger equations satisfied by the contour Green function are 

[dl + ml{x)]A{x, y) = -6^^\x, y) + f d^x'U{x, x')A{x', y) , (10) 

•J o 

[dl + mi{y)]A{x, y) = -6^'^ (x, y) + J^ rfVA(x, x')U{x' , y) , (11) 

where n(x, y) is the self-energy; the integration over Xq is performed on the contour 
and the function (5*^^)(x, ?/) is defined on the contour as 

5*^^^ {x — y) for Xq , I/O from the upper branch, 
5^^\x,y) = \ for xo , y^ from the different branches, 

—5'-^-' (x — y) for Xq , yo from the lower branch. 
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Let us split the self-energy into three parts as 

n(x,?/) =Us{x)5^'^\x,y) + Il>{x,y)e{xo,yo) + I\.<ix,y)Q{yo,xo) . 

As we shall see later, Us provides a dominant contribution to the mean-field while 
n< determines the collision terms in the transport equations. Therefore, we call Us 
the mean-field self-energy and n< the collisional self-energy. 

With the help of the retarded and advanced Green functions (||) and the retarded 
and advanced self-energies defined in an analogous way, the equations (p!0[) and ( |TT1) 
can be rewritten as 

[dl + ml{x) - Us{x)]A"<{x,y) 

= fd^x'[u^{x,x')A~{x',y) + U+{x,x')A^{x',y)], (12) 

[dl + ml - Us{y)]A^{x,y) 

= fd^x'[A^{x,x')U-{x',y) + A+{x,x')U^{x',y)], (13) 

where all time integrations run from — oo to +cxd. 

Let us also write down the equations satisfied by the functions A^ 

[dl + mlix) - Usix)]A^ix,y) 

= -(5(^)(x -y) + I rfVn±(a;, x') A±(x', y) , (14) 

[dl + ml{y) - Usiy)]A^ix,y) 

= -5'^^\x-y)+fd''x'A^{x,x')U^{x',y). (15) 

5 Towards Transport Equations 

The transport equations are derived under the assumption that the Green func- 
tions and the self-energies depend weakly on the sum of their arguments, and that 
they are significantly different from zero only when the difference of their arguments 
is close to zero. To express these properties it is convenient to define a new set of 
variables as 

A{X,u) = AiX + ^u,X-^u). 

For homogeneous systems, the dependence on X = (x + y)/2 drops out entirely due 
to the translational invariance and A{x,y) depends only on u = x — y. For weakly 
inhomogeneous, or quasihomogeneous systems, the Green functions and self-energies 
are assumed to vary slowly with X. We additionally assume that the Green functions 
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and self-energies are strongly peaked near u = 0. The effective mass m^,{x) is simply 
assumed to be weakly dependent on x. 

We will now convert the equations ( p!2| , |13|) into transport equations by implement- 
ing the above approximation and performing the Wigner transformation (^ for all 
Green functions and self-energies. This is done using the following set of translation 
rules which can be easily derived: 



d x'f{x,x')g{x',y) 
h{x)g{x,y) 

h{y)9{x,y) 

dUfix^y) 
d^f{x,y) 



fiX,p)giX,p) 



+ 



dfiX,p) dg{X,p) dfiX,p) dg{X,p) 



dp^ 
h{X)g{X,p)- 



dXt" dXt" dpf, 

idh{X)dg{X,p) 



h{X)g{X,p) + 



2 dXi" dp^ 
t dh{X)dg{X,p) 
2 dXt' dp^ 






Here X = {x + 2/)/2, d^ = -t^ and the functions f{x,y) and g{x,y) satisfy the 
assumptions discussed above. The function h{x) is assumed to be weakly dependent 



on X. 



Applying these translation rules to eqs. (|1^, [T^), we obtain 



1^2 _ ^pM^^ _ p2 ^ ml{X)-Us{X)-'-d,{ml{X)-Us{X))di:]A'<{X,p) 
= Il^{X,p)A-iX,p) + U+{X,p)A^{X,p) 
+ l{lf<{X,p), A-iX,p)} + ^{n+(X,p), A^(X,p)} , (16) 



^d' + zp^d,-p' + ml{X)-Us{X)+'-d,{ml{X)-IlsiX))di:]A'^{X,p) 
= A^<{X,p)U-iX,p) + A+{X,p)U^{X,p) 

+ '-{A^<ix,p), n'ix,p)} + ^{a+(x,p), n^(x,p)} , (17) 

where we have introduced the Poisson-like bracket defined as 

dC{X, p) dD{X, p) dC{X, p) dD{X, p) 



[C{X,p),D{X,p)] 



dp. 



^X^' 



^X^' 



dp. 



The kinetic theory deals only with averaged system characteristics. Thus, one usu- 
ally assumes that the system is homogeneous on a scale of the Compton wave length 



of the quasiparticles. In other words, the characteristic length of inhomogeneities 
is assumed to be much larger than the inverse mass of quasiparticles Therefore, we 

impose the condition 

1 



A<{X,p) 



> 



mt 



;d'A<{X,p) 



li 



which leads to the quasiparticle approximation. As discussed in the next section and 
in the Appendix, the requirement (0) renders the off-shell contributions to the Green 
functions A< negligible. Thus, we deal with the quasiparticles having on-mass-shell 
momenta. Unfortunately, the assumption (|1^) cannot be applied to massless particles 
and for this reason we have introduced the effective mass m^^. 

Let us now take the difference and the sum of eqs. (|TBp and (|TT), where the d^ 
terms have been neglected due to the quasiparticle approximation (18|). Then, one 
gets 



'p^d, + -d,{ml{X) 



Yis{X))d^]A^{X,p) 

^(n>(x,p)A<(x,p) - n<(x,p)A>(x,p; 
^{n+(x,p) + n+(x,p), A^(x,p)} , (19) 



-/ + ml{X)-Iis{X)\A^{X,p) 

= ^(n^(x,p)(A+(x,p) + A-(x,p)) + (n+(x,p) + n-(x,p))A^(x, 
+ '-{n>{x,p), A<(x,p)} - ^{n<(x,p), A>(x,p)} , 



p 



(20) 



where we have used the identity (0) applied to the Green functions and self-energies. 
One recognizes eq. (PD| ) as a transport equation while eq. (PDj ) as a so-called mass- 
shell equation. We will write down these equation in a more compact way. From the 
definition (IBt) one finds that 



A±(X,p) = ±i(A>(X,p) - A<(X,p) 



2ni J ijj — u)' 

The first term in the r.h.s is antihermitian while the second one is hermitian. Thus, 
we introduce 



ImA±(X,p) = ±i(A>(X,p) - A<(X, 



V. 



(22) 



ReA^{X,p) 



2m 



uj — UJ' 



(23) 



With the help of (^) and analogous formulas for 11^ the equations (|1^, ^0|) can 
be rewritten as 

[p'-ml{X) + U,{X) + Ren+(X,p), A^(X,p)} 

= z(n>(x,p)A<(x,p) - n<(x,p)A>(x,p)) 

- {n^(X,p), ReA+(X,p)}, (24) 



p2 -m2(X) + n5(X) + Ren+(X,p)]A^(X,p) = -n^(X,p)ReA+(X,p) 

{n>(x,p), A<(x,p)} + ^{n<(x,p), A>(x,p)} . (25) 



I 
4 



In the case of fields with the finite bare mass the gradient terms in the right-hand- 
sides of eqs. (p^ , p5l) are small [0, |^ and are usually neglected. When the bare fields 
are massless, as those studied here, there is no reason to neglect the gradient terms. 
The equation analogous to (|^) was earlier derived in 0, H- 

It appears useful to write down the transport and mass-shell equations satisfied 
by the retarded and advanced Green functions. Starting with eqs. (0,0) one finds 



{p2 - ml{X) + ns{X) + n±(X,p), A±(X,p)} = 



(26) 



p 



ml{X) + Tls{X) + n±(X,p) A±(X,p) = 1 . 



(27) 



We observe that the gradient terms drop out entirely in eq. (pTj). Nevertheless the 
equation holds within the first order of gradient expansion. Due to the absence of the 
gradients, eq. (pTf) can be immediately solved as 



(28) 



A±(X,p) 



p2-m2(x) + n,(x) + n±(x,p) 



One notices that A^ of the form ( pSf ) solves not only eq. ( [27| ) but eq. ( pG]) as well. 
Indeed, any function f of K satisfies the equation {K, f{K)} = 0. 

The real and imaginary parts of A''=(X, p), which are needed in our further con- 
siderations, are 



ReA±(X,p) 



p^ 



m2(X) + n,(X) + Ren+(X,p) 



(p2 - ml{X) + UsiX) + Ren+(X,p))' + (Imn+(X,p))' ' 



(29) 



ImA^(X,p) 



±Imn+(X,p) 



(p2 - ml{X) + Us{X) + Ren+(X,p))' + (Imn+(X,p))^ 



(30) 
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6 Free quasiparticles 



Before further analysis the equations obtained in the previous section we consider 
here a very important hmit which corresponds to the free quasiparticles. Specifically, 
we assume that II^ = 11^ = 0. Then, the equations (24, 25) read 



p'-d, + ^d,ml{X)d^]Al{X,p) = 



(31) 



[p'-ml]Al{X,p) = 0. (32) 

Although the quasiparticles are assumed to be free, the transport equation is of the 
Vlasov not of the free form. This is a simple consequence of the X— dependence of 
the effective mass. 

Due to eq. (|3^) Aq{X,p) is proportional to 6{p'^ — ml), and consequently free 
quasiparticles are always on mass-shell. If the quasiparticle approximation ( pTSD is 
not applied, the mass-shell equation gets the form 

'^d'-p' + ml{X)]AUx,p) = 0, 

and the off-shell contribution to the Green function Aq is nonzero. A detailed dis- 
cussion of the quasiparticle approximation is given in the Appendix. 

We also discuss the (anti-)chronological Green functions Ag in the limit of free 
quasiparticles. One easily finds their equations of motion as 



1 
2' 



p''d, + -d,mt{X)d;\A',{X,p) = 

1 . 



[p'-ml]A^,{X,p) 

For the antichronological function A", the right hand side of the mass-shell equation 
equals —1 instead of +1. The solution of these equation can be written as 



AS(X,p) 



p2 _ ^2 _|_ ^QH 



+ ei~po)A>{x,p) + eipo)A<ix,p) 



where Aq{X,p) is assumed to satisfy eqs. (|3T|, ^). It is worth mentioning that any 
function which depends on {X,p) through (p^ — ml) solves the Vlasov equation (|3T|). 
Aq{X,p) obeys the initial condition of the standard Feynman propagator. It also 
satisfies the relation (H). 

The antichronological Green function is 



AS(^,P) 



-1 



+ Q{po)A>{X,p) + Q{-po)A<{X,p) 



p2 _ ^2 _ ^g 

Knowing Aq and Aq one immediately gets the retarded and advanced functions 

1 



Ao(^,P) 



p2 _ ^2 -|- j^p^Q 



+ 



(33) 
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which obey the respective initial conditions. Confronting the expressions of A^ for 
free ( p3D and interacting quasiparticles (^), one finds that 

Imm(A-,rt = -Imn-(.Y,p) = {°: ["J P» > ° ' (34) 

in the hmit of free quasiparticles. 

It appears useful to express the Green functions Ag through the distribution 
function /o as 

Q{po)iA<{X,p) = e(po) 27r6{p^-ml)fo{X,p) = I^S{E, - po)fo{X,p) , (35) 



where Ep = wp^ + m^ . This equation should be treated as a definition of /q. 
Due to the relation (^ we have 

A<{X,p)=A>{X,-p), 

and consequently 

e(po)iA>(X, -p) = e(po) 2n6ip^ - mDUX, p) = ^5(i?p - Po)/o(X, p) . (36) 

In that way we express the positive energy part of Aq{X,p) and the negative 
energy part of Aq (X, p) through fo{X, p). We extend these expressions to the whole 
energy domain using the identity (|^). With the help of the explicit form of the 
retarded and advanced functions (^) we get the formula 

tA>{X,p) - iA<{X,p) = 27rSip' - ml){Q{j>^) - e(-po)) , (37) 

which is discussed in detail in the next section. 

Combining eqs. (^, ^6]) with (|37D one finds the desired expression of the Green 



functions Ag in terms of the distribution function /q. Namely, 

iA<iX,p) = ^SiEp - po)foiX,p) + ^S{Ep + po)ifoiX,-p) + 1) , (38) 

iA>iX,p) = ^5iEp-po)ifo{X,p) + 1) + ^S{Ep + po)fo{X,-p) . (39) 

-C/p rjp 

When the system is in thermodynamical equilibrium the distribution functions 
reads 

where (3^^ = u'^/T with m'^ being the hydrodynamical velocity and T the temperature. 
In the local equilibrium the two parameters are X— dependent. 
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7 Spectral function 

In this section we introduce one more function which appears useful in the anal- 
ysis of the interacting systems. The spectral function A is defined as 

A{x, y) = iA> {x, y) - i A< (x, y) . 

Thus, 

A{x,y) = mx),<l){y)]), 

where [0(x),0(?/)] denotes the field commutator. 
Due to the equal time commutation relations 

[0(t, x), 0(t, y)] = , [0(t, x), 0(t, y)] = -z<5(3) (x - y) , 

with the dot denoting the time derivative, the Wigner transformed spectral function 
satisfies the two identities 

I^JlA(X.p)=0. /^p„^(X,rt = l. (41) 

One also sees that (cf. eq. (p2D) 

A{X, p) = t2 ImA± (X, p) . (42) 

Finally, we observe that the identity (|^), which holds for the real fields, provides the 
relation 

A{X,p) = -A{X,-p). 

From the transport and mass-shell equations (|2^ , ^) one immediately finds the equa- 
tions for y4(X,p) as 

{p2-m^(X) + n5(X)+Ren+(X,p), A{X,p)} = 2{lmn+(X,p), ReA+(X,p)} , (43) 

[p^ - ml{X) + n5(X) + Ren+(X,p)] A(X,p) = 2Imn+(X,p) ReA+(X,p) . (44) 

Substituting ReA+(X, p) from eq. ( P^ ) into the algebraic equation (|3^) we find 
its solution as 

MX V) = 2lmU^X,p) 

^ '^^ {p2-ml{X) + Us{X) + ReU+{X,p)f + {ImU+{X,p)f' ^ ' 

Then, one easily shows that the function of the form (|45| ) solves eq. ( |4^ ) as well. In 



fact, the spectral function (^) can be found directly from eq. ([29| ) due to the relation 
(0)- 
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The spectral function of the free quasiparticles can be, obviously, obtained from 
]) but the limit should be taken with care. We firstly write the spectral function 
]) as 

AiX,p) 



p2 - ml{X) + Us{X) + ReU+{X,p) + ilmU+{X,p) 



p2 - ml{X) + Us{X) + ReU+{X,p) - iImU+{X,p) ' 



Then we take the limit 11 — > keeping in mind the condition (plf). Using the well- 
known identity 

I ■n+ =P-Ti'r^S{x) , 
X ± zO+ X 

we get the spectral function of noninteracting quasiparticles as 

Ao{X,p) = 2n6ip^ - ml){e{po) - 0(-po)) , (46) 



which, of course, coincides with 

Let us also consider a specific approximate form of the spectral function. If the 
condition 

p2 + ml{X) - UsiX) - Ren+(X,p) > |Imn+(X,p)| 

is satisfied, A as a function of po is close to zero everywhere except two narrow regions 
around po = ^E± {E± > 0) which solve the equations 

EliX, p) = p' + ml{X) - Us{X) - Ren+(X, po = ±E±, p) . 

Then, the function (03) can be approximated as 

1 r+(x,p) 



E+{X,p) (E+(X,p)-po)' + Tl{X,p) 

1 r_(x,p) 

E_(X,p) {E.{X,p)+pof + Tl{X,p) 
with 



(47) 



^ ^^ 2E±(X,p) ^ 



One easily checks that the spectral function of the form {^7\) satisfies the sum rules 

8 Perturbative expansion 



As discussed in e.g. R y, 11| the contour Green functions admit a perturbative 



expansion very similar to that known from vacuum field theory with essentially the 
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same Feynman rules. However, the time integrations do not run from — oo to +00, 
but along the contour shown in Fig. 1. The right turning point of the contour (tmax) 
must be above the largest time argument of the evaluated Green function. In prac- 
tice, to is shifted to — cxd and tmax to +cxd. The second difference is the appearance 
of tadpoles, i.e. loops formed by single lines, which give zero contribution in the 
vacuum case. A tadpole corresponds to a Green function with two equal space-time 
arguments. Since the Green function A{x,y) is not well defined for x = y we ascribe 
the function —iA^{x,x) to each tadpole. The rest of Feynman rules can be taken 
from the textbook of Bjorken and Drell ||2^ . 



In this section we consider the lowest-order contributions to the self-energies. It 
should be stressed that the Green functions, which are represented by the lines of the 
Feynman diagrams, correspond to those of free quasiparticles not of noninteracting 
fields. 

0^ model 

The lowest-order contribution to the self-energy which is associated with the 
graphs from Fig. 2 is 

U{x,y) = 6'-^\x,y){ml{x) - — A<(x,x)) , 

giving 

ns{x)=ml{x)-'-^A<{x,x), (49) 

and 

U>{x,y) = U<{x,y) = 0. 

Substituting A^ given by eq. (pS]), where A^ is expressed through the distribution 
function, into (|49|) one finds II^ as 



n.W = ,n;w-|/-i?^|2/„(x.p) + i|. 



As seen the integral is quadratically divergent in the vacuum limit when /o — * 0. This 
type of divergence, which appears due to the zero-mode fluctuations, is well known in 
the field theory. We remove it away by subtracting the vacuum value from II^. Thus 
one gets 

(2^p- 
We compute 11^ for the equilibrium system when the distribution function is given 
by eq. (|D) with u^ = (1,0,0,0). For T ^ m^ we get after elementary integration 



n.w = ™;w-|/-^/„(.,p). (50) 



the well-known result, see e.g. |I7[ , 



2, ^ 9T\x) 



Us{x) = mtix) - '-^ , (51) 

where the temperature is x— dependent in the case of local equilibrium. 
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model 



The lowest-order contribution to the self-energy corresponding to the graphs from 
Fig. 3a and 3b is 



Iia{x,y) = 5^'\x,y) 



ml--g J d x'Ao{x',x)A^{x',x') 



Locating the argument x on the upper branch of the contour, one finds 



Us{x) =ml- -/ / dV|A^(x',x) - A<(x',x) 



A<(a;',x') 



where the time integration runs from — oo to +cx3. Observing that A^ — A< 
we get 

Tlsix) 



A+ 



^l - 2^' / ^VA+(x',x)A<fx',x'^ 



Using the explicit form of Aq" ( p3D and expressing A"^ through the distribution 
function, one finds 



nAfx) 



g^ f d^x'd^p d?k 



-ip{x'—x) 



'^*^''' 2 7 (27r)4 {2n)^Ek p^-m^, + tpoO^ 



fo{x',k), (52) 



where as previously we have subtracted the vacuum contribution. 

In the case of global equilibrium, when fo{x, k) is independent of x, the integral 
from (|52|) can be computed as 

^2 " d^k 



^5 = ml + ^ I ^,,r__/Q(k). 



2ml J (2vr)3Efc 



In the limit T ^ m^, one finally finds 

n5 = 



fn^ H ^ • 

* 24 m2 



(53) 



(54) 



The graph from Fig. 3c corresponds to 



Ub{x,y) = --g'^Ao{x,y)Ao{y,x) 



and it gives 



U<ix,y) 



-yA<ix,y)A>iy,x). 



Substituting A< expressed through the distribution function to eq. 



(55) 
one finds 



d^k 



d^q 



(27r)^5(=^)(p + q-k) 



n (X,p) ^g J ^27r)32Efe (27r)32E, 

'6ipo-Ek-Eg)fo{X,k)fo{X,-q) 

+ 6ipo -Ek + E,) /o(X, k) (/o(X, q) + 1) 

+ 6{po + Ek- E,) (/o(X, -k + 1) /o(X, -q) 

+ 6{po + Eu + E,) (/o(X, -k) + 1) (/o(X, q) + 1) 



(56) 
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' 6{po -Ek- E,) (/o(X, k) + 1) (/o(X, -q) + 1) 

+ 5{po -Ek + Eg) ifoiX, k) + 1) /o(X, q) 

+ 6{po + Ek- Eg) /o(X, -k) (/o(X, -q) + 1) 

+ 6{po + Ek + Eg)fo{X, -k)/o(X, q)] . (57) 

Let us observe here that due to the trivial kinematical reasons H^ (X, p) = U"^ {X, p) = 
for on-mass-shell four-momenta p i.e. when p^ = ml. 

We further compute Il^{X,p) for p = 0. In the case of equilibrium when the 
distribution function is of the form (j^) and, as previously, u^ = (1, 0, 0, 0), one finds 
n^(X,p) for p = {±uj,0) with u; > as 

n>(x,cj,o) = n<(x,-c<;,o) 

e(^-2m,)Wl-— i- -2, (58) 



167r V ^ 



(exp(f)-l 



n<(x,cj,o) = n>(x,-Lj,o) 



-—Q{uj-2m,)\l ^. 59 

1671 V uj^ rexp(f)-l 



We also compute 

imn+(x,±cu,o) = ^(n>(x,±cu,o)-n<(x,±cj,o 



2i 



6(u; — 2m*) 



167r ' *' V 0^2 gxp (^) - 1 ■ 

One sees that the condition ( |5^ is indeed satisfied by the perturbative self-energy 
when g —>■ and m* is kept fixed. 



9 Quasiparticle mass 



The structure of eqs. (|2^ , [25| ) motivates the definition of the quasiparticle mass 
as a solution of the equation 

n5(X) + Ren+(X,po = m„p = 0) = 0. (60) 

The definition is not Lorentz invariant but statistical systems usually break such an 
invariance. In the case of global equilibrium there is, for example, a preferential 
reference frame related to the thermostat. 

Let us now look for the explicit expression of m^ in the lowest nontrivial order in 
g within the two models under consideration. 
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(j)^ model 

The only nonvanishing contribution to the self-energy in the lowest order of 
perturbative expansion conies from lis, which is given by eq. (|50|) . Therefore, the 
effective mass is a solution of the equation 

2/ A 9 f d^P 



■m^(x) = — - — -- — fn(x,p) . (61) 

*^ ' 2 J {2TxfEp^ '^' ^ ' 

One should keep in mind that Ep and /o depend on m,,. 

For the equilibrium system with T ^ m,, we immediately get the well known 
result: 

ml{x) = ^ . (62) 

One sees that the condition T ^ m^ is automatically satisfied in the perturbative 
limit where g^^ ^ 1. 

(j)^ model 

Now we have the nonzero contributions to the self-energy not only from II^ ( p2D 
but from 11 < (^, ^7^ as well. The real part of 11''' which enters eq. (|60| ) is given by 



the equation analogous to (pS]). Thus, 

T. ^+/.. N It./", n>(x,cj,o)-n<(x,Lj,o) 

Ren+ X, m*, = P / du ^ ' ' ^ k , , J _ 

2m J m^ — ui 

Using eqs. (B3, ISO) one finds 



Ren+(X,m„0) = --l-mj^ -^Jl ^""^ 



327r2 
where 



167r2 * J2m, uj^-ml^ uj^ exp (f^) - 1 

^ I{(3m,) , (63) 



^^> - A x2 - 1/4 e"- - 1 ■ 



Thus, in the case of global equilibrium with T ^ m^ the equation ( |60D reads 

24 ml 327r2 






One observes that 

,, , 1 r°° dx 1 
1(a) < - — = - . 



a Ji x"^ a 



Thus, I{[3m^) < T jm^ and consequently the absolute value of the third term in l.h.s 
of eq. (0) is much smaller then that of the second one in the limit of T ^ m^. 
Therefore, there is no real positive m* which satisfies eq. (p^). It means that the 
massive quasipartides do not emerge in the massless (p^ model. It is not surprising 
since the potential V{(j)) = y/P^ from the lagrangian (|I]) has no minimum, even a local 
one. So, we do not consider the 0^ model any more and concentrate on the 0^ model. 



Having the mass of quasipartides we can determine their dispersion relation. Since 
all self-energies except Ils{X) vanish in the lowest order of the perturbative expansion, 
the mass-shell equation ( pS]) coincides with that one of the free quasipartides (|3^). 
Thus, we have explicitly shown that in the first order of g, the 0^ model provides the 
system of free quasipartides described by the transport equation of the Vlasov form 
i.e. 

P^d^ + ^d,ml{X)di:]fo{X,p) = , (65) 

with the effective mass given by eq. (|6l|). 

10 Higher order self-energy 

We discuss here the g"^ contributions to the self-energy in the 0^ model, which 
are represented by graphs shown in Fig. 4. The bubble is again related to the effec- 
tive mass. The contributions corresponding to the diagrams 4a and 4b can be easily 
computed. However, they are pure real and the only effect of these contributions 
is a higher order modification of the effective mass. Thus, we do not explicitly cal- 
culate these diagrams but instead we analyse the contribution 4c which provides a 
qualitatively new effect. 

The graph from Fig. 4c gives the contour self-energy as 

^c{x,y) = —Ao{x,y)Ao{y,x)Ao{x,y) , 



and consequently 



„2 

^^ix,y) = ^A^{x,y)Al{y,x)A^{x,y) . 



Substituting Ag expressed through the distribution function /o as in eqs. (^) and 
(^) one finds the self-energy H< in the form 

U<(X,p) = t^ f , "^'^ , "^y , '^ll (27r)^^(^) (p + q - k - r) 
^ '^' 6 J 27rY2E, 27iY2E„ 27iY2Er.^ ' ^^ ^ ' 



6{po-Ek-E,-Er)ftfo'f^ 
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+ S{po -Ek-Eq + Er) f^ fo' iU' + 1) 
+ 6{po-E, + Eq-Er) f^ ifS + l) r, 
+ 5{p, -E, + Eq + Er) f^ ifS + l)(/o"- + 1 
+ 6{po + E,-Eq- Er) if,-' + 1) W fo 
E 



+ 6{po + E,-Eq + Er) iW + 1) fo' UV + 1) 

+ 5ipo + E, + Eq- Er) Uo' + 1) (/o + 1) /o 

+ 5(po + E, + Eq + Er) Uo' + 1) (/o + l)(/o"' + 1) 



(66) 



n>(x,p) 



(fk 



d?q 



d?i 



(27r)V3)(p + q-k-r) 



6 J (27r)32Efc {2Tif2Eq {2nY2E, 

S{po + Ek + Eq + Er)fo''fSfo' 
+ d{po + Ek + Eq- Er) U' fo (/o + 1) 
+ 8(po + E,-Eq + Er) U' ifo' + 1) W 
+ S{po + Eu-Eq- Er) fo" ifo' + l)(/o + 1) 
+ 5{po -Eu + Eq + Er) Uo + 1) /o fo'' 
+ 8{po -E, + Eq- Er) (/o^ + 1) f^ (/^ + 1) 
+ 5{po -E,-Eq + Er) {f'o + 1) Uo' + 1) /o"^' 
+ 5{po -E,-Eq- Er) Uo + 1) Uo' + l)(/o + 1) 



(67) 



with 



/o^^/o(X,k), /-'=^/o(X,-k). 

It is important to notice that in contrast to the similar expressions of the 0^ model 
i.e. eqs. ( p6|) and (^), the self-energies ( |67D and ( |66D are nonzero not only for the 
off-shell but for on-shell momentum p as well. However, the number of terms from 
eq. (|67|) or (|66D which contribute to n< is reduced when p^ = ml. Indeed, eqs. 
simplify in this case as 



e(po)n<(x,p) 



.r 



Sk 



d^q 



<Pi 



6 J {2'Kf2Ek {27T)^2Eg {2T:f2E, 

5{po ~E,-Eq + Er) f^ fo' ifo' + 1) 
+ S{po -Eu + E,- Er) f'o m + 1) /o 



(27r)^,5(3)(p + q-k-r) 





^{X^P) 


+ 5{j)o + Ek- 
.g"^ f cPk 


- Eg — Er 

(fq 

{2n)^2Eg 

d^q 
(27r)32E, 
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) ifo' + 1 
d^r 


' Jo Jo 
{2ti)^5^^\p + q - 

(27r)^5(=^)(p + q- 


-k- 
-k- 


r\ 


e(po)n= 


~ ^ 2 J {2n)^2Ek 
X (/o' + l)/oVo, 

.g^ r (Pk 


{2n)^2Er 
d^r 


- r) 
(68) 


~ \ J {2n)^2Ek 


{2'K)nEr 


^) 



Sipo -Ek-E, + Er) (/o' + 1) {fo' + 1) fo' 
+ 5{po -Ek + E,- Er) (/o' + 1) n (/^ + 1) 
+ 5{po + E,-E,~ Er) fo' ifo' + 1) (/o + 1) 



d^k 



d^q 



d^r 



■3_ 

^ 2 ] {2Tif2Ek (27r)32E, {2Tif2E, 

fo Wo + 1) (/o + 1) , 



{2T:f5'^^\p + q-k-r) 



(69) 



where po = y^* + P^- 

The self-energies n< provide, through the equations analogous to (|22D and (|23|) , 
the Ren"*" and Imll"*" which enter the transport (^) and mass-shell (^) equations. 
The imaginary part of 11+ is of particular interest. Due to the finite value of Imn+, 
the spectral function ( ^Sj) is no longer delta-like but it is of the Breit-Wigner shape. 
Thus, the quasiparticles are of finite life time. For the on-mass-shell momenta with 
Po > the imaginary part of 11+ equals 



Imn+(X,p) 



1 
2i 

4 



n>(x,p)-n<(x,p) 



d'^k 



d^q 



d^7 



X 



(27r)32Efe {27r)^2Eg {27t)^2E, 
/oM/o+l)(/o + l)-(/o + l)/o/o 



{2ny5^^\p + q-k-r) 



(70) 



This function was computed for the equilibrium distribution in ^7\, |T^] 



11 Interacting quasiparticles 

In this chapter we discuss the dispersion relation of the interacting quasiparticles 
and then define the respective distribution function. 

Having the self-energies calculated in g^ order we can determine the quasiparticle 
dispersion relation in this order. As previously, the quasiparticle mass is found as a 
solution of the equation (|60D. Therefore, the singular self-energy Ils{X) as well as 
Ren"'"(X, m*, 0) are included in ml{X). To avoid the double counting, the expression 



p 



is everywhere replaced by 



with 



ml{X) + Us{X) + ReU-'{X,p) 



p^-ml{X) + ReU+{X,p) 



Ren+(X,p) = +Ren+(X,p) -e(po)Ren+(X,m„0) -0(-po)Ren+(X,-m„O) 
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The dispersion relation is given by the equation 

p'-ml{X) + Ref[+{X,p) = 0., (71) 



but according to eq. (^51) , which determines the spectral function, the quasiparticles 
are of the finite width and the relation ([71| ) gives only the most probable quasiparticle 
four- momentum. We call the four-momenta, which satisfy the relation (0), as 'on- 
mass-sheir, however one should keep in mind that the meaning of this term differs 
for the finite and zero width quasiparticles. 

Eq. ([7l| ) can be easily solved if Rell"*" provides only a small correction to the free 
quasiparticle dispersion relation. Then, one finds the on-mass-shell momentum as 
p^ = i±E^,p) with 



E^ = \jml{X) + p2 + Ren+(X, ±^ml{X) + p^, p) . (72) 

The distribution function f{X,p) of the interacting quasiparticles is defined in a 
way analogous to eq. (|35| ) i.e. 

Q{po)iA<{X,p) = e(po) A{X,p) f{X,p) , (73) 

where A{X,p) is the spectral function (^). In contrast to the case of free quasiparti- 
cles, f{X,p) depends not on the three- vector p but on the four- vector p. Due to the 
identities 

A<{X,p) = A>(X, -p) , A{X,p) = iA>{X,p) - iA<{X,p) , 

we have 

iA>iX,p) = e(po) AiX,p) if{X,p) + 1) - 0(-po) A{X,p) f{X, -p) , (74) 

iA<{X,p) = e(po) A{X,p) fiX,p) - e(-po) A{X,p) {fix, -p) + 1) . (75) 



There is a very important property of A^ expressed in the from ( [f4t [f5|). Namely, if 
the Green functions A< satisfy the transport equation (p^ and the spectral function 
solves the equation (|4^) , the mass-shell equation of A<, i.e. eq. (|25|) , is satisfied 
automatically in the 0-th order of the gradient expansion. Let us derive this result. 

The transport and mass-shell equations ( p4t pSf ) with the gradient terms neglected 
read 

= n>(x,p)A<(x,p) - n<(x,p)A>(x,p) , 

"/ -m^(X) + Ren+(X,p)]A^(X,p) = -n^(X,p)ReA+(X,p) . 

Substituting A^ in the form ( |74| , [75| ) into the first equation and taking only the terms 
corresponding to po > 0, the equation is 



= A{x,p) 



n>(x,p) /(x,p) - n<(x,p) (/(x,p) 
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(76) 



Now we substitute A< given by ( |75D into the mass-shell equation and get 

p2 _ ml{X) + Ren+(X,p)l A{X,p) f{X,p) = ~U< {X , p)ReA+ {X , p) , (77) 



where po is assumed to be positive. Using the spectral function equation (|4^) , the 
equation ([77|) is manipulated to the form 



ReA+(a;,p) 



U>{X,p) f{X,p) - n<(X,p) (/(X,p) + 1 



. (78) 



One sees that if / solves eq. (^), it automatically satisfies eq. (|78D . Similar consid- 
erations can be easily repeated for po < and then for A^ with po > and po < 0. 

Let us observe that the quasiparticles studied in this paper are narrow. Indeed, the 
effective mass (|6T| ) is of order g^^'^ while the width of the Breit-Wigner distribution 



given by eq. (|48| ) is proportional to g in power at least 3/2. Thus, the width of 
quasiparticles is much smaller than their mass {g is obviously assumed to be small). 
Due to this property we will often refer to the case of zero-width quasiparticles or 
on-mass-shell momenta. 



12 Transport equation 

The distribution function / satisfies the transport equation which can be obtained 
from eq. (p4|) for A^ or A^. After using eq. (|43|) one finds 



A{X,p){p'-ml{X) + Ren+(X,p), /(X,p)} 

= tA{x,p) (n>(x,p) f{x,p) - n<(x,p) (/(x,p) + 1) 

+ z/(X,p){n>(X,p), ReA+(X,p)} 

- t{fiX,p) + 1) {n<(X,p), ReA+(X,p)} , (79) 

where po > 0. We have also used here the following property of the Poisson-like 
brackets: 

{A, BC} = {A, B}C + {A, C} B . 

Since eq. (^) is one of the main results of this paper we discuss it in detail. 

The left-hand-side of eq. (^) is a straightforward generalization of the drift term 
of the standard relativistic transport equation. Computing the Poisson-like bracket 
and imposing the mass-shell constraint ([71]) one finds the familiar structure 



^e{po){p'-ml{X) + Ren+(X,p), /(X,p 

d 



KVqI + vv)/(x,p) + vv{x) v,f{x,p) 



where 

V{X) = ml{X)-ReIl+{X,p) 
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and the velocity v equals dE^ /dp with the energy E^ given by eq. (|7^). 

Let us now analyse the right-hand-side of eq. (^). Since the quasiparticles of 
interest are narrow, we take into account only those terms contributing to the self- 
energies n^ which are nonzero for the on-mass-shell momenta. The other terms are 
negligibly small. Then, n< from the transport equation ( [75| ) are given by the formulas 
analogous to eqs. (|68|,|69D with / instead of /q. Consequently, 



Q{po)U<{X,p) 



g^ f d^kA+ d^qA+ d^rAt 



2 J (27r)4 (27r)4 

{r + 1) f r , 



{2nf6^^\p + q-k- 



(80) 



e{po)u>{x,p) 



g^ f d^kAt d^qA+ d^rAt 



2 J (27r)4 (27r)4 (27r)4 

X /M/' + 1) (r + 1) , 



(27r)^(5(^)(p + 



^1) 



where 



At^Q{ko)A{X,k). 
One sees that in the limit of zero-width quasiparticles 



d^kAl 

72^ 



d'^k 



{2'Kf2Ek 



The first term in r.h.s of the transport equation ([79|) is very similar to the standard 
collision term of the relativistic transport equation pTI . Indeed, 



z(n>(x,p)/(x,p) - n<(x,p)(/(x,p) + i) 

g^ f d^kAl dSA+ d^rAl^ 



(82) 



{2ixf5^^\p + q~k~ 



X 



2 J (27r)4 (27r)4 (27r)4 
(/p + 1) (/<? + 1) f f _ fP fi ^f + 1) ^f + 1) 



The last two terms from r.h.s of eq. ([79|), analogous to those found long time ago 
in [§, 1^, are absent in the usual transport equation. We are going to show that in 



the local equilibrium, when the collision term (|82D vanish, we reproduce the standard 
collisionless equation if the four-momentum is on mass-shell. 
As well known ^1 



the standard collision term, which emerges from (|8^ ) when 
the quasiparticle width tends to zero (cf. eq. (|46|)), vanishes for the local equilibrium 
distribution function of the form 
collision term ( 



q). Following [21 



one easily shows that the 
also vanishes for the distribution function (BOl) with the particle 



momentum p no longer constrained by the mass-shell condition. 

One observes that in the local equilibrium the coUisional self-energies can be writ- 
ten as 

U>{X,p) = 2tlmU+{X,p) {r\X,p) + 1) , 
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and 

n<(X,p) = 22lmn+(X,j9) r\x,p) . 

The transport equation ([79|) then simphfies to 

A{X,p){p'-ml{X) + Ren+(X,p), r^(X,p)} 

= 2Imn+(X,p) {r{X,p), ReA+(X,p)} . 

Using eqs. (^) and (^5]) one manipulates this equation to the form 

Imn+(X,p){p2 - ml{X) + ReIl+{X,p),r{X,p)} (83) 

= (/ _ ml{X) + Ren+(X,p)){lmn+(X,p), r^(X,p)} . 

As seen, the term on the r.h.s drops down for the on-mass-shell momenta and then 
we reproduce the usual Vlasov equation i.e. 

{p' - ml{X) + Ren+(X,p), r{X,p)} = . (84) 

The role of the two unusual terms from r.h.s of the transport equation ( ffOD beyond 
the local equilibrium is rather unclear and needs further studies. 

13 Summary and concluding remarks 

We have discussed in this paper the nonequilibrium features of the massless 
fields. The derivation of the kinetic equation in such a case faces serious difficulties 
because there is no natural length scale over which the system inhomogeneities can 
be integrated over. As known the transport theory deals with the quantities averaged 
over an elementary phase-space cell of the minimal size given by the particle Compton 
wavelength. 

The fields with the zero bare mass usually gain an effective mass due to the self- 
interaction. Therefore, we have introduced the auxiliary mass term in the lagrangian 
and then, the transport theory has been derived in a way very similar to the earlier 
studied 0, §] case of massive fields. However, due to the position dependence of the 
effective mass, the limit of the noninteracting quasiparticles corresponds to the Vlasov 
rather than the free particle case. The smallness of the effective mass has also forced 
us to take into account some extra gradient terms which are usually neglected in the 
transport equation. 

We have considered in detail the 0^ and 0^ models which appear to be very 
different from each other. In the 0^ model the effective mass is generated in the lowest 
nontrivial order of the perturbative expansion. In contrast, the massive quasiparticles 
do not emerge in the (p^ model and most probably there is no transport limit of this 
model which, as well known, is anyhow ill defined. 
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Within the 0^ model we have derived the transport equation for the finite width 
quasiparticles. The distribution function has been defined in such a way that the 
mass-shell constraint is automatically satisfied (in the gradient zeroth order). We 
have found except the mean-field and collision terms the specific ones which are 
absent in the standard transport equation. However, in the case of local equilibrium 
we have been able to reproduce the usual coUisionless equation if the four-mementum 
is on the mass shell. 

The massless fields play a crucial role in the gauge theories such as Quantum 
Chromodynamics. We believe that the methods developed in this study will be useful 
in the discussion of the transport theory of quarks and gluons. Till know only the 
mean-field limit of such a theory is well understood p2|, ^ . In spite of some efforts 



24, 25 1 , the systematic derivation of the collision terms is still missing. The problem 



of the off-mass-shell propagation, which on phenomenological level has been proven 



to play a very important role ||2^, is also unsettled. 
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Appendix 

We discuss here the quasiparticle approximation for the system of noninteracting 
fields. To simplify the discussion the bare mass m is assumed to be nonzero or 
equivalently m^, is treated as a constant. The transport equation and the mass-shell 
constraint read 

p-dAliX,p) = 0, (85) 

[jd^-p^ + m^]A^{X,p)=0. (86) 

These equations, which directly follow from the field equation of motion (Q) with 
m = m^,, are exact in the case of the massive free fields - the gradient expansion is 
not needed to derive them. 

The mass-shell constraint (|86|) shows that the function Ag {X, p) is indeed nonzero 
for the off-shell momenta i.e. Aq (X, p) 7^ for p"^ 7^ m"^. This result looks surprising 
if one keeps in mind that the field, which solves the equation of motion (|^), is in a 
sense on mass-shell. The field is the sum of the plane waves 



d^k 



[X 



f fe-'^" a(k) + e''^" a*(k)) , (87) 
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where k = (cu, k) with uj = -y/k^ + m^. Thus, k"^ = rn?. Substituting the field ( p?] ) into 
the A<(X, p) definition, one finds that the off-shell contribution to Ag (X, p) conies 
from the interference of the positive and negative energy parts present in eq. (|87|) . 
Let us consider when such a contribution can be neglected. 

One easily shows that the transport equation (|85D is solved by the function which 
depends on the four-position X = (t, x) only through x — vt i.e. 

A|(X,p) = F(x-vt,p). 
The quasiparticle condition (|T^) applied to the function F reads 



F(x — vt,p) 



> 



m^ 



,ViVj 



d^F{^-vt,p) 



■'^3^ 



(9(x — vt)id(x — vt)j 



If this condition is satisfied for every x at a given moment of time, say to, it is satisfied 
at any time. In other words, if the initial condition at to is sufficiently homogeneous 
that the quasiparticle approximation can be applied, then this approximation is ap- 
plicable at any time - the system remains homogeneous. 

The question arises whether Ag, which simultaneously solves the transport ( pSj ) 
and mass-shell ( p6D equations, can satisfy the quasiparticle condition. We introduce 
the Fourier transformed function Aq{Q,p) defined as 



AliQ^p)"^' Jd'Xe^^-''AUx,p) 



(88) 



The equations corresponding to (^) and (|86D , respectively, read 

p-qaUq,p) = o, 



[--Q^-p^ + m']Al{Q,p) = 0. 



They are both solved by 



tA<{Q,p) = 6{p-Q) 6{--Q'-p' + m') A{Q,p) 



(89) 



with A{Q,p) controlled by the initial condition. Since iAo(X, p) is real, A{Q,p) has 
the property 

A{Q,p) = A*{-Q,p). (90) 

The solution of the equations (|85|) and (|86|) satisfies the quasiparticle condition 
f|l8|) when 



A{Q,p) » %A{Q,p) 
or equivalently A{Q,p) ^ only for Q^ ^ rv? . 



(91) 
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It is instructive to consider the explicit solution of eqs. ( ]85| ) and ( ]86| ) in 1 + 1 
dimensions. Using ( pQ]) we get 



tA<{X,p) 



<PQ 



-iQ-X 



(27r)2 
0( V) + 0(P' - "^' 



5(p-Q)5(-Tg'-p' + ^')^(Q,p) 



pz 



(27r)^|p^| \ p^ — m? 



-iQX 



A(g,p) + e^«^A(-g,p) , 



(92) 



where Q denotes the two- vector 



Q.2|p„|,?^(^,l 

V P ^Po 

Keeping in mind the property (^) the solution ( p2D can be rewritten as 

zAUx,p) = [e(-/) + ©(/ - m')] {f{p) sin(gX) + g{p) cos(QX)) , (93) 

where f{p) and g{p) are the real functions of p determined by the initial condition. 
The quasiparticle condition ([ISD is satisfied by (|93D if 



One also sees that A^{X,p) ~ (5(p^ — m^) only for Q = 0. In other words, the 
function Ag {X, p) is strictly zero for the off-mass-shell momenta when the system is 
exactly homogeneous. If we are interested in the weakly nonhomogeneous systems, 
the functions are nonzero for p^ > rn?. Equivalently, if p^ = m? then p^ > rn^ but 
not p^ < m?. 

The properties of the function Ag in 1 -|- 1 dimension can be trivially generalized 
to 3 + 1 case showing the limitations of the quasiparticle approximation. 
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Figure Captions 

Fig. 1. The contour along the time axis for an evaluation of the operator expectation 
values. 

Fig. 2. The lowest-order diagrams of the self-energy in the 0*^ model. The bubble 
in a) denotes the additional interaction due to the effective mass. 

Fig. 3. The lowest-order diagrams of the self-energy in the 0^ model. The bubble 
in a) denotes the additional interaction due to the effective mass. 

Fig. 4. The second-order diagrams of the self-energy in the (p^ model. 
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